Impact of variance components on reliability of absolute quantification using digital PCR

Background

Digital polymerase chain reaction (dPCR) is an increasingly popular technology for detecting and quantifying target nucleic acids. Its advertised strength is high precision absolute quantification without needing reference curves. The standard data analytic approach follows a seemingly straightforward theoretical framework but ignores sources of variation in the data generating process. These stem from both technical and biological factors, where we distinguish features that are 1) hard-wired in the equipment, 2) user-dependent and 3) provided by manufacturers but may be adapted by the user. The impact of the corresponding variance components on the accuracy and precision of target concentration estimators presented in the literature is studied through simulation.

Results

We reveal how system-specific technical factors influence accuracy as well as precision of concentration estimates. We find that a well-chosen sample dilution level and modifiable settings such as the fluorescence cut-off for target copy detection have a substantial impact on reliability and can be adapted to the sample analysed in ways that matter. User-dependent technical variation, including pipette inaccuracy and specific sources of sample heterogeneity, leads to a steep increase in uncertainty of estimated concentrations. Users can discover this through replicate experiments and derived variance estimation. Finally, the detection performance can be improved by optimizing the fluorescence intensity cut point as suboptimal thresholds reduce the accuracy of concentration estimates considerably.

Conclusions

Like any other technology, dPCR is subject to variation induced by natural perturbations, systematic settings as well as user-dependent protocols. Corresponding uncertainty may be controlled with an adapted experimental design. Our findings point to modifiable key sources of uncertainty that form an important starting point for the development of guidelines on dPCR design and data analysis with correct precision bounds. Besides clever choices of sample dilution levels, experiment-specific tuning of machine settings can greatly improve results. Well-chosen data-driven fluorescence intensity thresholds in particular result in major improvements in target presence detection. We call on manufacturers to provide sufficiently detailed output data that allows users to maximize the potential of the method in their setting and obtain high precision and accuracy for their experiments.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2105-15-283) contains supplementary material, which is available to authorized users.     

Relationships between Sediment Microbial Communities and Pollutants in Two California Salt Marshes

Salt marshes are important ecosystems whose plant and microbial communities can alter terrestrially derived pollutants prior to coastal water discharge. However, knowledge regarding relationships between anthropogenic pollutant levels and salt marsh microbial communities is limited, and salt marshes on the West Coast of the United States are rarely examined. In this study, we investigated the relationships between microbial community composition and 24 pollutants (20 metals and 4 organics) in two California salt marshes. Multivariate ordination techniques were used to assess how bacterial community composition, as determined by terminal restriction fragment length polymorphism and phospholipid fatty acid analyses, was related to pollution. Sea urchin embryo toxicity measurements and plant tissue metabolite profiles were considered two other biometrics of pollution. Spatial effects were strongly manifested across marshes and across channel elevations within marshes. Utilizing partial canonical correspondence analysis, an ordination technique new to microbial ecology, we found that several metals were strongly associated with microbial community composition after accounting for spatial effects. The major patterns in plant metabolite profiles were consistent with patterns across microbial community profiles, but sea urchin embryo assays, which are commonly used to evaluate ecological toxicity, had no identifiable relationships with pollution. Whereas salt marshes are generally dynamic and complex habitats, microbial communities in these marshes appear to be relatively sensitive indicators of toxic pollutants.     

Diversity, composition, and geographical distribution of microbial communities in California salt marsh sediments

The Pacific Estuarine Ecosystem Indicators Research Consortium seeks to develop bioindicators of toxicant-induced stress and bioavailability for wetland biota. Within this framework, the effects of environmental and pollutant variables on microbial communities were studied at different spatial scales over a 2-year period. Six salt marshes along the California coastline were characterized using phospholipid fatty acid (PLFA) analysis and terminal restriction fragment length polymorphism (TRFLP) analysis. Additionally, 27 metals, six currently used pesticides, total polychlorinated biphenyls and polycyclic aromatic hydrocarbons, chlordanes, nonachlors, dichlorodiphenyldichloroethane, and dichlorodiphenyldichloroethylene were analyzed. Sampling was performed over large (between salt marshes), medium (stations within a marsh), and small (different channel depths) spatial scales. Regression and ordination analysis suggested that the spatial variation in microbial communities exceeded the variation attributable to pollutants. PLFA analysis and TRFLP canonical correspondence analysis (CCA) explained 74 and 43% of the variation, respectively, and both methods attributed 34% of the variation to tidal cycles, marsh, year, and latitude. After accounting for spatial variation using partial CCA, we found that metals had a greater effect on microbial community composition than organic pollutants had. Organic carbon and nitrogen contents were positively correlated with PLFA biomass, whereas total metal concentrations were positively correlated with biomass and diversity. Higher concentrations of heavy metals were negatively correlated with branched PLFAs and positively correlated with methyl- and cyclo-substituted PLFAs. The strong relationships observed between pollutant concentrations and some of the microbial indicators indicated the potential for using microbial community analyses in assessments of the ecosystem health of salt marshes.     

Soil Microbial Community Associated with an Invasive Grass Differentially Impacts Native Plant Performance

This study is one of the first to show that invasive plant-induced changes in the soil microbial community can negatively impact native plant performance. This greenhouse experiment tested whether soil microbial communities specific to the rhizospheres of an invasive grass (Aegilops triuncialis) and two native plants (Lasthenia californica and Plantago erecta) affected invasive and/or native plant performance. Each of these species were grown in separate pots for 2 months to prime the soils with plant-specific rhizosphere microbial communities. Each plant species was then planted in native- and invasive-primed soil, and effects on plant performance were monitored. At 5 months, differences in microbial biomarker fatty acids between invaded and native soils mirrored previous differences found in field-collected soil. L. californica performance was significantly reduced when grown in invaded soil compared to native soil (flowering date was delayed, aboveground biomass decreased, specific root length increased, and root mass ratio increased). In contrast, P. erecta and A. triuncialis performance were unaffected when grown in invaded vs native soil. These results suggest that in some cases, invasion-induced changes in the soil microbial community may contribute to a positive feedback loop, leading to the increased dominance of invasive species in an ecosystem.     

Detritus, trophic dynamics and biodiversity

Traditional approaches to the study of food webs emphasize the transfer of local primary productivity in the form of living plant organic matter across trophic levels. However, dead organic matter, or detritus, a common feature of most ecosystems plays a frequently overlooked role as a dynamic heterogeneous resource and habitat for many species. We develop an integrative framework for understanding the impact of detritus that emphasizes the ontogeny and heterogeneity of detritus and the various ways that explicit inclusion of detrital dynamics alters generalizations about the structure and functioning of food webs. Through its influences on food web composition and dynamics, detritus often increases system stability and persistence, having substantial effects on trophic structure and biodiversity. Inclusion of detrital heterogeneity in models of food web dynamics is an essential new direction for ecological research.     

Application of Real-Time PCR To Study Effects of Ammonium on Population Size of Ammonia-Oxidizing Bacteria in Soil

Ammonium oxidation by autotrophic ammonia-oxidizing bacteria (AOB) is a key process in agricultural and natural ecosystems and has a large global impact. In the past, the ecology and physiology of AOB were not well understood because these organisms are notoriously difficult to culture. Recent applications of molecular techniques have advanced our knowledge of AOB, but the necessity of using PCR-based techniques has made quantitative measurements difficult. A quantitative real-time PCR assay targeting part of the ammonia-monooxygenase gene (amoA) was developed to estimate AOB population size in soil. This assay has a detection limit of 1.3 × 10⁵ cells/g of dry soil. The effect of the ammonium concentration on AOB population density was measured in soil microcosms by applying 0, 1.5, or 7.5 mM ammonium sulfate. AOB population size and ammonium and nitrate concentrations were monitored for 28 days after (NH₄)₂SO₄ application. AOB populations in amended treatments increased from an initial density of approximately 4 × 10⁶ cells/g of dry soil to peak values (day 7) of 35 × 10⁶ and 66 × 10⁶ cells/g of dry soil in the 1.5 and 7.5 mM treatments, respectively. The population size of total bacteria (quantified by real-time PCR with a universal bacterial probe) remained between 0.7 × 10⁹ and 2.2 × 10⁹ cells/g of soil, regardless of the ammonia concentration. A fertilization experiment was conducted in a tomato field plot to test whether the changes in AOB density observed in microcosms could also be detected in the field. AOB population size increased from 8.9 × 10⁶ to 38.0 × 10⁶ cells/g of soil by day 39. Generation times were 28 and 52 h in the 1.5 and 7.5 mM treatments, respectively, in the microcosm experiment and 373 h in the ammonium treatment in the field study. Estimated oxidation rates per cell ranged initially from 0.5 to 25.0 fmol of NH₄⁺ h⁻¹ cell⁻¹ and decreased with time in both microcosms and the field. Growth yields were 5.6 × 10⁶, 17.5 × 10⁶, and 1.7 × 10⁶ cells/mol of NH₄⁺ in the 1.5 and 7.5 mM microcosm treatments and the field study, respectively. In a second field experiment, AOB population size was significantly greater in annually fertilized versus unfertilized soil, even though the last ammonium application occurred 8 months prior to measurement, suggesting a long-term effect of ammonium fertilization on AOB population size.     

The glucose transport system of the hyperthermophilic anaerobic bacterium Thermotoga neapolitana

The glucose transport system of the extremely thermophilic anaerobic bacterium Thermotoga neapolitana was studied with the nonmetabolizable glucose analog 2-deoxy-D-glucose (2-DOG). T. neapolitana accumulated 2-DOG against a concentration gradient in an intracellular free sugar pool that was exchangeable with external source of energy, such as pyruvate, and was inhibited by arsenate and gramicidin D. There was no phosphoenolpyruvate-dependent phosphorylation of glucose, 2-DOG, or fructose by cell extracts or toluene-treated cells, indicating the absence of a phosphoenolpyruvate:sugar phosphotransferase system. These data indicate that D-glucose is taken up by T. neapolitana via an active transport system that is energized by an ion gradient generated by ATP, derived from substrate-level phosphorylation.     

Ustilago maydis Mating Hyphae Orient Their Growth toward Pheromone Sources

Snetselaar, K. M., Bolker, M., and Kahmann, R. 1996. Ustilago maydis mating hyphae orient their growth toward pheromone sources. Fungal Genetics and Biology 20, 299-312. When small drops of Ustilago maydis sporidia were placed 100-200 μm apart on agar surfaces and covered with paraffin oil, sporidia from one drop formed thin hyphae that grew in a zig-zag fashion toward the other drop if it contained sporidia making the appropriate pheromone. For example, a2b2 mating hyphae grew toward a1b1 and a1b2 mating hyphae, and the filaments eventually fused tip to tip. Time-lapse photography indicated that the mating hyphae can rapidly change orientation in response to nearby compatible sporidia. When exposed to pheromone produced by cells in an adjacent drop, haploid sporidia with the a2 allele began elongating before sporidia with the a1 allele. Sporidia without functional pheromone genes responded to pheromone although they did not induce a response, and sporidia without pheromone receptors induced formation of mating hyphae although they did not form mating hyphae. Diploid sporidia heterozygous at b but not at a formed straight, rigid, aerial filaments when exposed to pheromone produced by the appropriate haploid sporidia. Again, the a2a2b1b2 strain formed filaments more quickly than the a1a1b1b2 strain. Taken together, these results suggest that the a2 pheromone diffuses less readily or is degraded more quickly than the a1 pheromone.